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We analyze the Sivers asymmetry in light-cone gauge. The average transverse momentum of the 
quark distribution is related to the correlation between the quark distribution and the transverse 
component of the gauge field at x~ ± oo. We then use finiteness conditions for the light-cone 
Hamiltonian to relate the transverse gauge field at = ±oo to the color density integrated over 
x~ . This result allows us to relate the average transverse momentum of the active quark to color 
charge correlations in the transverse plane. 
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O ' I- INTRODUCTION 

(N : 

^ , Many high energy inclusive hadron processes show surprisingly large transverse polarizations or single-spin asym- 
, O ■ metries The most prominent example is inclusive hyperon production, but similar phenomena are observed in 
' many other reactions as well. For example, in the inclusive photo-production of pions on a transversely (relative to the 
I ■ photon momentum) polarized nucleon target, a left-right asymmetry of the produced pions is observed. Theoretically, 
two mechanisms (which are not exclusive!) have been identified as a potential source of the asymmetry: the Sivers 
and the Collins mechanisms. In the Collins mechanism the asymmetry arises when a transversely polarized quark 
Cn| fragments into pions with a left-right asymmetry. In contrast, in the Sivers mechanism 3] the asymmetry results 
^ from an intrinsic transverse momentum asymmetry of the quarks in the target nucleon. At first, such an intrinsic 
' transverse momentum asymmetry was expected to vanish due to time-reversal invariance of the strong interaction. 
1^ \ However, more recently it was realized that, even at high energies, the final state interactions (FSI) of the struck quark 

■ play an important role for the single-spin asymmetry Formally the FSI can be taken into account by introducing 
' an appropriate Wilson line gauge link along the trajectory of the ejected quark The gauge invariantly defined 
1 transverse momentum distributions with a gauge link along the light-cone to infinity are no longer required to vanish 

■ to to time-reversal invariance and a nonzero Sivers asymmetry is possible. 
However, while the above reasoning explains the existence of the Sivers asymmetry it leaves many questions unan- 

O • swered, for example what sign and what magnitude should one expect for the asymmetry, i.e. is it just an obscure 
I ' small effect or is it large? If there is a large asymmetry, what does any information about the asymmetry teach 
us about the structure of the nucleon? In fact, it is also possible that there is no simple connection between the 
asymmetry and ground state properties of the nucleon since the asymmetry hinges on the inclusion of FSI. In this 
paper an attempt will be made to make a step towards answering these questions. 

The paper is organized as follows: in Section 2, we review the definitions of gauge invariant transverse momentum 
distributions and the role of the gauge field at — ±oo in the light-cone gauge. In Section 3, we use finiteness 
^ ■ constraints for light-cone Hamiltonians to derive operator constraints that allow one to relate the gauge fields at 
[ x~ = ±oo to degrees of freedom at finite x~ . In Sections 4 (QED) and 5 (QCD) we use these operator constraints to 
relate the average transverse momentum to charge (color charge) correlations in the transverse plane. 
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II. INITIAL (FINAL) STATE INTERACTIONS AND TRANSVERSE SPIN ASYMMETRIES 

Ref. explains how final state interactions (FSI) and initial state interactions (ISI) allow the existence of T-odd 
parton distribution functions. Formally the FSI (ISI) can be incorporated into kj^ dependent parton distribution 
functions (PDFs) by introducing a gauge string from each quark field operator to infinity ^ 

g(x,k^,S^) = |^?^-i!Zie--P+?'-+»k,-y-(p|5-(y-,y^)^+[y-,y^;00-,yJ [^^ (1) 

We use light-front (LF) coordinates, which are defined as: — (y+, y^, yj,), with j/^ — {y^±y^)/^/2. In all correlation 
functions, y^ = and we therefore do not explicitly show the dependence. The path ordered Wilson-line operator 
from the point y to infinity is defined as 

[oo^,y±;y^,y±] = Pcxp(^-ig J dz^A+{y+,z^,y±)^ . (2) 
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FIG. 1: Illustration of the gauge links in Eq. J^. 



The specific clioice of patli in Eq. |^ reflects the FSI (ISI) of the active quark in an eikonal approximation. The 
complex phase in Eq. is reversed under time-reversal and therefore T-odd PDFs may exist |5|, which is why a 
nonzero Sivers asymmetry |^ is possible. 

Naively, the single spin asymmetry seems to be absent in light-cone gauge A'^ = 0, since the Wilson lines in Eq. 
are in the x~ direction and therefore J dz^A^ = 0. Without the phase factor any single spin asymmetry vanishes 
due to time reversal invariance. 

This apparent puzzle has been resolved in Ref. |H, where it has been emphasized that a truly gauge invariant 
definition for unintegrated parton densities requires closing the gauge link at = oo, i.e. a fully gauge invariant 
version of Eq. JQ) reads 

g(x,k^,s^) - 1 ^y-fy^e-^^p^y-+^^-y- (3) 
X {p \q{y^,y±h^ [y^:y±;oo^,y±] [oo",y±, 00^,0^] [oo",o_l;0",o^] q{Q)\p) . 

In all commonly used gauges, except the light-cone gauge, the gauge link at x~ = cxd is not expected to contribute 
to the matrix element, since the gauge fields are expected to fall off rapidly enough at oo. However, this is not true 
in light-cone gauge and therefore it has been suggested in Ref. 6] that, in light-cone gauge, the entire single-spin 
asymmetry arises from the phase due to the gauge link at = oo 

q{x,k^,s^) = I dyJy± ^-^.p^y-+^^.^■y^ {p\q{y-,y^) [oo' , , ^ , 0±] 7+g(0)|p) ■ (^+=0) (4) 

This implies for the average transverse momentum for quark flavor q 

J d^k^qix,k^,s^)k^ = -gj ^e-'-P^y- {p\q{y-,0^)-f+A^i^-,0^)qiO)\p) (A+ = 0) (5) 

A similar result holds for the unintegrated parton density relevant for the Drell-Yan process, where the initial state 
interaction provides the asymmetry 

qpastix,k^,s^) = y'^^-i?Zie--P^^"+*-y-(p|g(y-,y^)[-oo-,y^,-oo-,0j7+g(0)|p) (A+ = 0) (6) 

and 

J d'k^qpastix,k^,s^)k^^-g J ^e-"P+^"(p|g-(y-,O^)7+Ax(-(^-,O^)g(0)|p) (A+ = 0) (7) 

(Here and in the rest of this paper we will work in light-cone gauge and therefore no longer emphasize explicitly that 
A~^ = 0). The fact that these asymmetries hinge on the value of the transverse gauge field at x^ — ±oo makes the 
evaluation of these matrix elements rather tricky. Only a careful regularization prescription for the — singularity 
of the gauge field propagator is capable to generate the complex phase that is necessary for a non-vanishing SSA. 
Therefore, the question arises whether knowledge of the light-cone wave function for the nucleon would (in principle) 
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sufficient information to calculate the SSA nonperturbatively, since the abovementioned prescription is perturbatively 
defined. 

One of the main purposes of this paper is to render Eqs. ^ and Q into a more useful form. For this purpose we 
first use the fact that the asymmetric part of the unintegrated parton densities relevant for semi-inclusive DIS and 
Drell-Yan are equal and opposite ^ 



Therefore 



k±q{x) = / (fk^q{x,k±,s±)k± 



(8) 



(9) 



J d'^kj_q{x,kj^,s^)kj^ - J d^k^qpast{x,kj^,Sj^)kj 



dy 
47r 



For later applications, we also rewrite Eq. @ in the more familiar color component notation 



q{y-,0^h+^q{0)al{0^: 



2 J An 

where Aq are the Gell-Mann matrices and 

<(b_L) = Ai,(cx)",b_L) - Ai^(-cx)",bj^) 



(10) 



(11) 



Eqs. (PI and (|10|l still involve the gauge field at x — ±oo. In the following section, we will derive an operator 
relation that relates A_L(±(X),y^) to degrees of freedom at — oo < x^ < oo. 



III. FINITENESS CONDITIONS 



In a light-cone formulation of QCD (x^ is time), not all components of the quark and gluon field are dynamical, 
since they satisfy certain constraint equations (Lagrangean contains no x'^ derivative for these degrees of freedom). 
In particular, A~ as well as the "bad" component of the quark field qhad = 57^7~9 are such constrained degrees of 
freedom in QCD. Upon eliminating both and q^ad using their constraint equations, one arrives at the canonical 
light-cone Hamiltonian for QCD in light-cone gauge — 0, which governs the x~^ evolution 



P- 



dx^d^yn I F 



-F. 



-tr 



where 



^12 



diA2-d2Ai+ig [^1,^2] 
i [dx + idy + ig (A^ + iAy)] u- 



(12) 



(13) 



i [dx - idy + ig {A^ 



J = d^d^A, + J 



iAy)] u+ 



V2 

m 

T2 



u+ 



d-d,A, + ig [A„d^A^] + g (^q+ql + 



(14) 



In the above expression q± form respectively the positive and negative chirality components of the 'good' component 
qgood = of the quark field. A summation over quark flavors is implicit. 

The requirement that acting on a hadron state is free of infrared x~ = ±00 divergences implies that each of 
the four terms in Eq. I|12() are free of such divergences. This has a number of consequences. In Ref. |9|, it was 
investigated what conditions on the hadron Fock state result from the requirement that such divergences are absent 
in the first two terms. In this work we will focus on the third term. The condition 







j dx Tr 


-J 1 J 



— finite 



(15) 
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implies that [lo| 

/•OO 

dx^ J{x^ = Vx_L (16) 



This (weak) condition, which should hold as a condition on all physical states, forms one of the crucial ingredients of 
this investigation. Since the first term in J is a total derivative, Eq. (|16() implies (now expressed in terms of the more 
familiar cartesian components i — 1,2, which are implicitly summed when they appear in pairs) 



d^a^{x±) = —p{x^) = — J dx J(x ,xj 



(17) 



where 

q;*(x_l) EE ^'(cx)",x_l) - A'(-oo",x±) 
J(x-,Xi) = ig[A\d-A']+g(q+ql + q-q'_). (18) 

The physical meaning of p(xj^) is the total charge (quarks plus gluons) along a line with fixed x^. 
For later use, we also express the above results using color components (instead of matrix notation) 

9*q;^(x_l) = -paix±) = - J dx^ Ja{x^ ,x±) (19) 

J<,(x-,Xi) = -gfabcAld^Al + gY.qj+Y'^, (20) 

where Aq are the Gell-Mann matrices. In QED, the analogous conditions read and 

9*a*(x_L) = -p{ii^) = - J dx-J{x-,XA^) (21) 

J(x-,Xi) = Y,e,qj+q. (22) 
q 

An additional condition arises from the requirement that the tr (^^12)^ -term in the Hamiltonian is convergent at 
x^ = ±00: the field _L strength tensor itself must vanish at x^ — ±00 

.Fi2(±c»-,Xi) = 0, (23) 

i.e. {±oo^ ,x±) must be pure gauge 

AJ(+(X3-,Xi) = --y+(xi)9^"v:|(xi) (24) 
.9 

Ai{-oo-,K±) = --l/_(xi)o»^V'i(xi). 
9 

This allows us to gauge transform A^ {—(X)~ ,x±) to zero while preserving yl+ = 0. Since in this gauge A^(+oo~, x^) = 
is still pure gauge, i.e. A^ {oo~ = —^U^x±)d^U{xj_) with U{k±) = vl{x±)V+{xj_) we thus conclude that 

a'(x_L) = A'(oo-,x_l) - A'{-(X)-,x^) (25) 

must be (in this gauge) of the form 

a*(xi) = --C/^(xi)9't/(xi). (26) 
9 

Eq. H26(l together with Eq. 1)17(1 thus determine a^{x±) uniquely (up to some trivial constants). 

The above results have a number of applications as we will discuss below. First of all, Eqs. H19I21|I alert us again that 
in light-cone gauge one must not assume a vanishing of the gauge fields at x^ = ±00. However, the most important 
application of Eqs. H19I21() lies in the fact that it allows us to reexpress a^{x±) = A'(oo~,xj_) — A'^{—oo~,x±) in 
terms of other degrees of freedom. The interesting aspect about this observations is the fact that a'^{x±) also appears 
in the correlation function (j^ for the average transverse momentum. In the rest of this paper we will discuss the 
implication of this fundamental result. 
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IV. AVERAGE ± MOMENTUM IN THE ABELIAN CASE 

Before proceeding to QCD, we will first discuss an abclian theory, where the nucleon contains quarks of different 
flavor q with charges respectively as well as photons. For simplicity, we will consider here only the Sivers assymmetry 
averaged over all x 

(k^q) ^ J dx J d'^k^q{x,k^,s±^)k^ ^ J dxk^q{x). (27) 

Considering only the x-averaged asymmetry not only helps keep the resulting expressions simpler, but may also help 
to cancel possible divergences from endpoint singularities |l2j . 

In QED, the requirement that a*(x^) is pure gauge can be rewritten as 

a'(xj.) = -aV(xi), (28) 

where '/>(xj^) is some scalar function, which can be determined by solving the 2-dimensional Poisson equation. This 
yields 

aXx.)=-/^-^^p(y.). (29) 

What we have found is an operator condition (|29|l that determines the transverse gauge field at x~ = zLoo in terms 
of the charge density at —oo < x~ < oo. This condition needs to be satisfied for states in order to have an infrared 
finite energy. Upon inserting this result into the expression for {k±q) we find 

(k±,) = / ^-^{p\mi^mpiy^)\p) (30) 

W J |y_L| 

where 

p{y±) = E / dy-q'iy)i+q'{y) (3i) 

q' 

is the charge density (integrated over x^) from all flavors (actually, for symmetry reasons, q' = q does not contribute 
in Eq. 

Eq. H3l)|l is gauge invariant and provides a regularized expression for the Sivers asymmetry that depends on the light- 
cone wave function of the target only. Although this should be clear from our derivation, it should be emphasized that 
the final state interactions are included in Eq. Ij^UI) , but the gauge fields that give rise to the final state interactions have 
been reexpressed in terms of charge density correlations inside the target. The relevant density-density correlations are 
correlations in the transverse plane of the charge density integrated along x~ . The kernel in the correlation function 
is the Lorentz-boosted Coulomb force integrated along x~ as well. Such a result should not be surprising since this is 
simply the Coulomb force from the spectators acting on the escaping quark. In fact, for the specific example of the 
scalar diquark model this result was already obtained in Ref. [T^ . 

As a corollary, we should also note that if one sums the mean transverse momentum over all quark flavors one gets 
zero 



= 0. (32) 



This result follows from the fact that the integration kernel in Eq. 129|l is odd under x^ ^ y^. 

To summarize this section, what we have found is that the mean transverse momentum for quarks of flavor q can 
be related to the correlations between quarks of flavor q and all other quarks in the transverse plane. This is not 
surprising since the final state interaction is the Lorentz boosted Coulomb interaction and the correlations describe 
the Coulomb field from the spectators acting on the active quark. 

V. AVERAGE _L MOMENTUM IN QCD 

In QCD, the condition that the gauge field at ±oo is pure gauge is nonlinear, which prevents us from writing down 
closed form solutions to the finiteness conditions. This can be seen as follows. If one writes 

C/(xi) = e-*"'^"'^^)^ (33) 
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then, to lowest order in (j)a (see also Appendix B) one finds the QED-like condition 

A(/.,(xx) = -Pa(x±), (34) 

yielding 

|x_L - y_L| 



ali^^) = -aV, ^ - / ^Pa{y±). (35) 



Of course, there are nonabelian corrections to Eq. (|34|l and therefore, unlike in QED, Eq. I|35|l is not an exact solution 
to the finiteness conditions. However, since we were unable to find an exact operator solution, we will proceed using 
Eq. (123). 

Upon inserting H29() into Eq. H10() one obtains the nonabelian version of Eq. H30|l 



9(0)7+y9(0)p.(y^) 



(36) 



The physical interpretation of this result is that the average transverse momentum of quarks of flavor q can be 
related to correlations on the transverse plane. The specific correlations that appear in Eq. (|36|l reflect a Coulomb 
interaction between the active quark and the spectators. Here a Coulomb interaction appears because we have solved 
the finiteness constraints in QCD only to first order. Eq. (|36|l is thus equivalent to treating the FSI in lowest order in 
perturbation theory [l. ITllllsL IT^ . Note also that the resulting correlation fmictions are very similar to the correlation 
functions that have been used to describe the small- a; gluon distributions in nuclei . 

To first order, what we also find is that the average transverse momentum due to the FSI of all constituents (quarks 
+ gluons) added together vanishes for symmetry reasons. It is not clear if this happens beyond lowest order. 

Eq. (|36|l may be useful for several reasons. While the original expression for the Sivers asymmetry involved a gauge 
link, which made a parton model interpretation difficult, Eq. (??) does have an immediate parton model interpretation 
in terms of color-flavor correlations in the transverse plane. This may be useful in correlating experimental data with 
our understanding of the nucleon structure. Another use of Eq. 136|l is that it can be directly calculated from the 
light-cone wave functions of the nucleon. Of course, we need to keep in mind that (unlike the QED case) Eq. H3t)|) 
is only approximation, but we still believe that this result provides a step towards linking the Sivers asymmetry with 
other features of hadron structure. Finally, we would like to emphasize that Eq. (|36|) suggests interesting connections 
between the distribution of partons in impact parameter (x^ ) and the sign of the transverse SSA 

PI- 

For example, 

in a simple quark model, such as the bag model JT], the color part of the matrix element in Eq. (|35|l would be 
negative (attraction). If the transverse distribution of j+ is transversely shifted relative to the spectators (3 then 
the resulting average transverse momentum has the opposite sign of the sign of the transverse distortion in impact 
parameter space. 



VI. SUMMARY 



We have studied the average transverse momentum of gauge invariant Quark distributions for a transversely po- 
larized target in light-cone gauge. The Wilson line is along the light-cone to infinity to incorporate the Final state 
interactions in semi-inclusive DIS. In light-cone gauge, the Wilson-line phase factor receives its only nonzero con- 
tribution from the gauge field at — ±oo. In a naive Fock space expansion the gauge field at x~ = ±00 is 
usually implicitly set to zero, thus making a correct treatment of single-spin asymmetries rather difficult (except in 
perturbation theory, where one can carefully regularize the fields at x~ = ±00 "by hand" . 

We have also studied conditions for the infrared (x~ = ±00) convergence of the light-cone Hamiltonian for gauge 
theories and derived operator conditions that need to be satisfied in order for the light-cone to be free of infrared 
divergences arising from otherwise ill-defined operators . This operator condition relates the transverse component 

of the gauge field A^(±oo^,x^) to the color density pa{x±) integrated over all x~ 

Fortuitously, the same kind of operators that governs the average transverse momentum in gauge invariant quark 
distributions also appears in the finiteness conditions for the light-cone Hamiltonian. We are thus able to eliminate 
A^(±oo~,xj_) in the average transverse momentum in favor of other, less infrared singular, degrees of freedom. In 
QED we can solve the operator condition arising from finiteness conditions exactly and we are able to express the 
average transverse momentum in terms of charge density correlations in the transverse plane. In QCD we were only 
able to solve the operator condition to first order in the color charge density and there we find a similar result as 
in QED, namely that the average transverse momentum can be related to transverse correlations between the active 
quark and the spectators. 
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Single spin asymmetries do not have a simple parton model (or light-cone Fock space) interpretation. The main 
significance of our results is that we have found relations that allow to relate the average transverse momentum to 
operators that do have a parton interpretation (in QED exactly, in QCD approximately). One immediate application 
of these results is that it allows to evaluate the average transverse momentum of the quarks directly from the nucleon 
wave function in light-cone quark models. 

Several extensions of this work are conceivable. First it would be desirable to derive an exact solution (at least in 
terms of an expansion) for the finiteness conditions in QCD, so that one can study the effects of higher order terms 
that we have omitted. Secondly, it would be interesting to see if one can translate the results from the work into 
lattice language (Euclidean as well as transverse lattice) with the goal of being able to compute the average transverse 
momentum nonperturbatively within these frameworks. 

This work was supported by the Department of Energy (DE-FG03-96ER40965). I would like to thank J. Collins 
and X. Ji for useful discussions and comments. 

APPENDIX A: ALTERNATIVE DERIVATION FO THE FINITENESS CONDITION 

In this appendix we would like to present an alternative derivation of Eq. (|20(l that does not require a light-cone 
Hamiltonian. We start from the + component of the QCD equations of motion in light-cone gauge A'^ = 

D^Ti:+ ^ d-T-+ + d^T:+ - 9fabcKr+ - .7+ (Al) 

where 

T-+ = d-At - d+A- - gfabcA^At = -d^A' (A2) 
= d^A+ - d+A], - gfabcAlAt - -d.Al 

and the fermion color density is given 

j+^gY^q^+^q (A3) 
Q 

is the fermion color density. Written out in components the Field equations (jAl|) thus read 

- d'_A- - d^d^A\ - gfab.A',r+ = ,7+ (A4) 

Integrating Eq. ljA4|) over while making use of the condition that the field strength tensor at ~ ±oo 

vanishes, thus yields 



9M: 



/OO 
dx- [j+ - gfabcA'.d.A',] (A5) 
-OO 

dx" jQ(a;",xj_) = -/9a(x±) = /3(xi). 



APPENDIX B: HIGHER ORDER COLOR CORRELATIONS 

We want to solve the finiteness condition 

a^<(x^) - -p.(x^) (Bi) 

subject to the constraint that a\{ii.i_) is pure gauge, i.e. 

a'(xj.) = -^-U^d'U. (B2) 
9 

In order to satisfy the second condition, one can make the ansatz 

C/(x_l) = cxp (-ig(j)a{^i)^Y) ■ (^2) 



Upon inserting this ansatz into Eq. IjASp one finds 



a^(x^) = -aV(x^) - I [0(xx),9V(xi)] + 1^ [0(x^), [0(xx),9V(x±)]] + ... (B4) 
i.e. in component notation 

2 

a^(xi) = -5Va(x±) + pabcM^±WM^±) - Y^fabcM^±)fcdeM^±)d'de{^±) + ... (B5) 

and therefore 

9*<(xx) = -d'd'cba{^^) + PabcM^±)d'd'M^±) (B6) 
-J^fabcfcde [0h(Xi)<?^d(x_L)9'(9Ve(x±) + 9 Vh (x^ )<?id (x_L (x± )] + ... 

= -Pa(x±) 

One may attempt to solve Eq. Hij6|l by making a formal expansion in powers or pa(xj^) 

</),(x^) = ,^(°)(x^)+0W(x^) + ..., (B7) 

yielding 

5^aV(")(xx) = Pa(xi) (B8) 

fwr(xj9^av(°)(x^)^|. 



5^av(i)(x^) = f/,fc,</.i")(xj_)9*av(°)(x^) = |/„.fc,0i")(xi)p,(x^). 



e.t.c. 

Eq. (|35l) is then obtained by keeping only the lowest order term in Eqs. ljB5|) and 1|B8I 



[1] 
[2] 
[3] 
[4] 

[5] 
[6] 
[7] 
[8] 
[9 

[lo; 

[11 
[12 

[13 
[14 
[15 
[16 
[17 



P.J. Mulders, lectures at PRAHA 2001, hep-ph/0112225; D. Boer and J. Qiu, Phys. Rev. D 65, 034008 (2002). 
J.C. Collins, Phys. Lett. B 396, 161 (1993). 
D.W. Sivers, Phys. Rev. D 43, 261 (1991). 

S.J. Brodksy, D.S. Hwang, and I. Schmidt, Phys. Lett. B 530, 99 (2002); S.J. Brodsky et al., Phys. Rev. D65, 114025 
(2002). 

J.C. Collins, Phys. Lett. B 536, 43 (2002). 

X. Ji and F. Yuan, Phys. Lett. B 543, 66 (2002); A. Belitsky X. Ji, and F. Yuan, Nucl. Phys. B 656, 165 (2003). 
D. Boer, P.J. Mulders, and F. Pijlman, hep-ph/0303034. 

M. Burkardt, Adv. Nucl. Phys. 23, 1 (1996); S.J. Brodsky, H.C. Pauli and S.S. Pinsky, Phys Rept. 301, 299 (1998). 
F. Antonuccio, S.J. Brodsky, and S. Dalley, Phys. Lett. B412, 104 (1997). 

W.M. Zhang, in "Light-Front Quantization and Non-Perturbative QCD", Eds. J. Vary and F. Wolz, (IITAP, 1997), p. 141; 

M. Burkardt, in "New Directions in Quantum Chromodynamics" , AIP Conf. Proc. 494, p. 239 (1999); hep-th/9908195. 

F. Yuan, Phys. Lett. B575, 45 (2003). 

J. Collins, Acta Phys. Polon. B34, 3103 (2003). 

M. Burkardt and D.S. Hwang, hep-ph/0309072. 

A. Bacchetta, A. Schafer, and J.-J. Yang, hep-ph0309246. 

M. Burkardt, Phys. Rev. D 66, 114005 (2002); M. Burkardt, hep-ph/0302144. 

L. McLerran and R. Venugopalan, Phys. Rev. D 59, 094002 (1999). 

M. Burkardt, Int. J. Mod. Phys. A 18, 173 (2003); X. Ji, hep-ph/0304037. 



